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Abstract.  -We  consider^he  stability  of  difference  schemes  for  the  solution  of  the  initial  boundary 
value  problem  for  equation 

- -  u,  =  (A(x,  t)  ux),  +  B{x,  l)ux  +  C(x,  t)u  +  f(x ,  t), 

where  u,  A,  B,  C  and  /  are  complex  valued  functions.  Using  energy  methods,  we  establish  the 
stability  of  a  general  two  level  scheme  which  includes  Euler’s  method,  Crank-Nicolson’s  method 
and  the  backward  Euler  method.  If  the  coefficient  is  purely  imaginary,  the  explicit  Euler’s 

method  is  unconditionally  unstable,  For  this  case, -we- propose  a  new  scheme  with  appropriately 
chosen  artificial  dissipation,  which-we  proveKto  be  conditionally  stable. 
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1.  Introduction 

Finite  difference  methods  for  the  solution  of  equations  of  Schrodinger  type  have  been  studied  by 
many  authors,  and  extensively  applied  to  solve  practical  problems  in  many  disciplines  [1,  5,  6,  7,  8, 
9,  10,  11].  Since  conventional  explicit  schemes  are  unstable  [2, 4],  implicit  schemes  are  usually  used, 
especially  the  Crank-Nicolson  scheme.  In  [7],  D.F.  Griffiths  et  al  discussed  a  predictor-corrector 
scheme.  In  [8,  12]  the  existence  and  convergence  of  solution  for  this  difference  equation  have  been 
studied.  Since  the  solution  of  Schrodinger  equation  possesses  conservation  laws,  the  schemes  which 
satisfy  discrete  conservation  laws  have  also  been  investigated  extensively  [6,  7], 

We  consider  the  equation 

u(  =  {A[x,t)ux)x  +  B(x,t)ux  +  C[x,t)  u  +  /(x,t),  (1.1) 

in  the  domain  <?r(0  <  x  <  l,  0  <  t  <  T),  where  u(x,t),  A(x,t),  B(x,t),  C(x,t)  and  f{x,t) 
are  complex  functions,  and  R eA(x,t)  >0  and  |A(x,f)|  ^  0.  This  kind  of  equation  arises  in 
plasma  physics  and  acoustics  [1,  11].  Clearly,  (1.1)  involves  both  equations  of  Schrodinger  type  and 
parabolic  equations.  In  this  paper,  we  consider  the  initial-boundary  value  problem  for  (1.1),  with 
the  conditions  : 


«(«»o  =  2(x) 

(1.2) 

[x=0  =  U|*»J  =  0, 

(1-3) 

where  u(x)  is  a  complex  function. 

In  Section  2,  we  analyse  the  stability  of  a  general  two  level  difference  scheme  for  (1.1),  (1.2) 
and  (1-3),  which  is  a  direct  generalization  of  a  well-known  scheme  for  parabolic  equations,  and 
includes  the  explicit  Euler  scheme  and  the  implicit  Crank-Nicolson  scheme  as  special  members. 
Some  new  results  are  obtained  from  which  we  can  see  the  relationship  between  the  Schrodinger 
equation  and  the  parabolic  equation. 

If  R eA(x,t)  e  0,  the  above  mentioned  explicit  Euler  scheme  is  unstable.  Several  interesting 
stable  explicit  schemes  have  been  presented  in  [3]  for  the  case  of  the  simplified  equation  u/  =  tu**, 
and  applied  to  some  underwater  acoustics  problem  in  [2].  In  these  schemes,  Euler’s  method  is 
stabilized  by  appropriately  chosen  artificial  dissipation.  In  Section  4,  we  extend  these  results  to 
equation  (1.1). 

We  let  h  denote  the  spatial  mesh  size  and  divide  the  finite  interval  [0,1]  into  the  mesh  intervals 
by  the  points  x;  =  jh  ( j  =  0,1,...,  J),  where  Jh  —  l.  We  let  kn  denote  the  size  of  the  time 


step  at  the  n-th  step.  For  convenience,  we  shall  denote  u(xj,tn)  by  u".  We  also  use  the  following 
difference  operators  : 


. ,  u1+1  -  u? 

\  ~  Ut>  ~  kn  • 


D+uj  = 


_  Mj+l  -  Uj 


D.uj  =  2tza= i, 


Duj  = 


_  Uj+ 1/2  ~  — 1/2 


Duj  = 


Uj+ 1  -  ttj-l 
2h 


~(D+  +  D.)  tty. 


Hence,  we  have 

^(“i+i  -  2u;  +  u,_i)  =  D+D-uj  =  D2uj. 

[^y+i/2(«y+i  -  «y)  ~  A>-i/2(“>  ~  «/-x)]  =  D(ADuj). 

Finally,  for  any  function  <f>,  we  use  ^n+“  to  denote  a<f>n+l  +  (1  —  a)0",  for  0  <  o  <  1. 
Next,  we  give  our  definition  of  stability.  First  we  define  the  inner  product  for  u  and  v  : 


J- 1 

( u,v )  =  ^  ujtijh 
>« i 

where  v  denotes  the  complex  conjugate  of  v  and  the  norm  for  u  : 


INI  =  \/(u,u)  . 


Definition  1.1.  We  call  a  scheme  stable  if  the  solution  u"  satisfies  : 

H«"ii<ciii«oii+cJf:ii/,iit, 

1-0 

where  C\  and  Ci  are  constants  which  are  independent  of  n  and  h. 


2.  A  General  Two  Level  Scheme. 

We  are  going  to  consider  the  following  scheme  for  (1.1),  (1.2)  and  (1.3)  : 


u*?+1  —  u?  1 

J-1—L  -  - «r* )  -  -  *p?>i 

/u"+a  -  u”+a\ 

_  gn+<*  I  J  _  £"+<* un+a  _  rn+a 


j  *  1,2,...,J-  1 


(2.1) 


We  can  also  write  (2.1)  in  this  form  : 


u?+1  -  DAn+aDu^+a  -  B"+a  Du”+a  -  C?+°  u"+a  =  fV 
*)  j  )  j  j  )  J  i 


n+or 


(2.4) 


Clearly,  (2.1)  or  (2.4)  is  explicit,  implicit  and  the  Crank-Nicolson  scheme  when  a  =  0,  a  /  0  and 
a  =  5,  respectively. 

We  will  see  that  this  scheme  is  unstable  if  Re  A( x,  t)  =  0  and  a  =  0.  In  this  case,  we  will  give 
another  conditionally  stable  explicit  scheme  in  Section  3. 

To  facilitate  the  analysis,  we  first  transform  (2.1)  by  a  change  of  variable  to  eliminate  the  first 
order  term. 


Lemma  2.1.  Suppose  A  €  C3,  B  6  C2,  C  €  C1  and  Re  A  >  0,  |A|  >  ao  >  0,  then  we  can  choose 
a  function  <f>,  such  that  0  <  Mo  <  <  A/j,  j  «  0, 1,*  •  •,  J  where  Mo  and  Mi  are  constants 

and  under  the  transformation 

u"+a  =  ^+0v^+0  (2.5) 

(2.JJ  becomes  : 


„»+i  _ 


_  Jj  -  ^(«r"  -  v-r,«)] 

-  -  Hj+aVj+a  -  -  Q”+“r"+1-°  =  F”+° 


(2.6) 


where  Gj,  Hj,  Kj,  Qj,  Fj  are  bounded  in  Qr- 
Proof.  See  Appendix  A. 


We  now  state  our  main  stability  result  for  (2.6). 

Theorem  2.1.  Suppose  the  conditions  in  Lemma  2.1  are  satisfied.  If  £  <  o  <  1,  then  scheme  (2.1) 
is  stable. 

Proof.  In  order  to  establish  the  stability  of  (2.6),  we  are  going  to  estimate  ||vn+,||.  We  multiply 
(2.6)  by  0"+o  to  obtain  : 


(L(t>n+°),v"+“)  =  (F"+0,t>n+0)  . 


(2.7) 


Then  we  estimate  every  term  in  (2.7).  For  the  first  term,  we  have  : 

Re(vJ*+1,vn+°)  =  a  Re  (t?J?+1,  v"+1)  +  (1  -  a)Re(Vp+1,  vn) 

-eEiTP-^Ewi* 

-  fRe  E1'"0,"'  +  E1-"*1*" 

= 2ibz  nr  i*  -  e  ki2)  -  e  i«r‘  -  *"i2 

=  -rIM'll2  -  Ml2)  -  iL^£iMK"+,ll2- 

For  the  second  term,  we  have 

Re  (DAn+a Dv ,  t>n+“)  =  -Re  (An+° D-vn+a ,D_vn+“) 

=  —  ||\/Re  An+aD-vn+a  ||2. 

For  the  rest  of  the  terms  on  the  left  hand  side  of  (2.7),  we  have 

jgm+a  | 

icr  <r 1 «r  i  s  L^(i»r  12 + i«r  i2> 

|Fn+av«+a€n+«|  <  LX_i(|^+«|2  +  |f>»+“|2) 

wx ■©* >n  s  ^f^d-'^n2 + wrn 
wr^*-«ri s  i<?n  [<»  -  «)i<r  i + °wi]  [°i»r  i + <>  -  «)i«?i] 

S^O-T1!2 +  !<?!’)• 


So  we  have 


|Re  [I2crj-iv?-ivi+a + n  ^+arr + n  Kj+iv?tivj+a + 

<  Mi  [||»"+“||2  +  ||v"~||2  +  |K+"|I2  +  i||v"+,|l2  +  fll»“ll2 

=  3A/,||tin+ft||2  +  ^||v"+1||2  +  ^||«”||2, 

where  M2  is  a  constant  which  is  an  upper  bound  of  |C?|,  |£T|,  |A'|  and  |Q|. 

For  the  right  hand  side  of  (2.7),  we  get 

|/2e  (F”+°,  vn+0)|  <  ^||v"+®||2  +  ^||F"+“||2  . 
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Now  we  substitute  expressions  (2.8),  (2.9),  (2.10)  and  (2.11)  into  (2.7)  and  obtain 


-^-(!|u"+1!l1  -  Ml1)  -  +  l|VJteA«*i>-.*+*||> 

<  3M,|K+"||!  +  ^||»"+1||I  +  ^Ml1  +  i||v"+“||!  +  i||F"+»||! 


where 


<Af3(||vn+1||2  +  |kl2)  +  |||F"+“||2, 


M3  =  6M2  +  “g-  +  1  • 


(2.12) 


When  5  <  a  <  1,  from  expression  (2.12)  we  have 


rHllvn+1ll2  -  IKII2)  <  M3(ibn+1it2  +  IKII2)  +  dlfm+all5 


It  is  easy  to  see  that  when  k„  <  we  have 


K+'||!  <  ±±|^ll»T  +  2Mlf,”ll1 
<{l  +  8*„M3)||v"f  +  2);„||r-||2. 


According  to  Lemma  4.1  in  Appendix  B,  we  obtain 


ll*n+Ml2  <  2c8A'*‘"+‘  IM12  + 


Av .ii,.' i i; ly  Cocl^s 
;  Aw .  H-^i/or 


(2.13) 


From  (2.5)  and  0  <  M0  <  \4>j\  <  M\  we  obtain 


||tin+,||I  <  [iM1  +  Dl/T*'  ■ 


(2.14) 


Now  we  turn  to  discuss  the  case  0  <  ot  <  5. 


Theorem  2.2.  If0<a<  ^  and 


o/i  (i-i  f  I  m>x  \An+°\7  1  1 

2(1  2a)  \  r  jfcjzr*  +  t, 


where  rj  is  an  arbitrary  positive  constant  then  the  scheme  (2.1)  is  stable. 


(2.15) 


«  y 


Proof.  First  let 


^=^iA-r/2iK++r-ri 

1  i  * n+Q  |  |.,n+or  .,n+or| 


rv = -p\A'x/2\\vra  - 
P3;  =  |67^IK-l°i 
p<j=\Hr\K+a\ 

i*-wr*ihn 

From  (2.6)  we  obtain  immediately 

i«p;'iJ<  (e«,)  • 

Next  we  construct  a  quadratic  form  in  Pij  (/  =  1,  •  •  • ,  7)  by  : 


(2.16) 


Pj  = 


Ijn+a  |2  1;  '  l^n 

'A/+l/2l  r*i 


i.jn+a  7  /  7  \‘ 

E  j?  -  a  -  E  ^ 

2-1/21  ;«3  \/«i  / 


(2.17) 


where  jj  is  an  arbitrary  positive  constant.  According  to  Lemma  4.2  (in  Appendix  C),  Pj  is  nonneg¬ 
ative  if 

kn  £ 


(1  -  2a) 

'  UST/sP  , 

kjw3 

+  I 

[***•*}£» 

h>R*A;±?/3 

+  n 

(2.18) 


which  is  true  because  of  (2.15).  Making  this  assumption,  we  have  from  (2.16)  : 


,  „  (h2ReA"+°2\  ,  (h2ReA"+°\  ,  2  , 

(i  -  i  <  )  pv  +  (  +  5,SP' 


Summing  up  these  expressions  from  j  —  1  to  j  —  J  —  1,  we  obtain 

(1  -  2a)kn\\v?+ltf  <  2||v/jfe ’a^D-v"+°\\2  +  5T?Af|(3||t;"+*||2  +  ||t>w+1-“||2)  +  59||F"+*||2. 
<  2||\/j?e  An+aD~ v"+“||2  +  4CM/22(|K+1||2  +  ||t>"||2)  +  5^||F"+°||2 


(2.19; 


Substituting  (2.19)  into  (2.12),  we  have 


ir(l|t,"+1||2  "  IKI!2)  <  ^4(||vn+I||2  +  Kll2)  +  Ms||F"+“||2, 
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where  A/4  =  2  A/3  +  40//A/22  and  A/5  =  1  +  5  r}.  It  is  easy  to  see  that  when  kn  <  5^-,  then 

<  (1  +  4fc„A/4)j|t,"||2  +  2A/5*n|lF"+l2  . 


According  to  Lemma  4.1,  we  obtain 


||wB+1||2  <  e4^'"+1(||n°||2  +  2 A/5  £  ll^+“ll2^) 


(2.20) 


and  from  (2.5)  and  0  <  Mq  <  |<6|  <  A/,  we  obtain 


1/2  ^ 

!«n+1||2  <  ^e‘,A/4'*+l(ll«0||2  +  2A/5£||F'+“||2*/) 


(2.21) 


which  implies  that  the  scheme  (2.1)  is  stable. 


3.  A  stable  explicit  scheme  for  the  case  Re  A  =  0. 

In  this  section,  we  assume  Re  A  =  0  and  let  a(x,  t)  denote  the  imaginary  part  of  A.  We  see 
from  the  previous  section  that  the  scheme  (2.1)  with  a  =  0  is  unstable.  This  is  unfortunate  because 
in  many  applications  an  explicit  scheme  is  desirable  because  they  tend  to  be  easier  to  implement, 
especially  in  a  vector  or  parallel  computing  environment.  In  [2],  we  construct  a  stable  explicit 
scheme  for  the  simple  equation  u<  =  iuzz  by  adding  appropriately  chosen  artificial  dissipative 
terms  to  (2.1).  Here  we  consider  an  extension  of  this  scheme  for  the  more  general  equation  (1.1). 
We  construct  the  following  scheme  : 

«7+1  -  v?  ~  1  (fa?  -  |o?|) 

- 1 - taJ - U - 4 - K+1  2v;  +  . 

*n  n  *  (3.1) 

-  -  °iu?  =  fh  J  =  1 . 1,  • 


where 


VJ  ~  rjK+i  ~  2«"  +  J  =  1, 


i’o  =  vj  —  0. 

We  will  prove  that  this  scheme  is  conditionally  stable.  For  convenience,  we  assume 


fly  >  0,  j  —  0, 1, . . . ,  J. 


(.  .  da  d2a  .  dB  .  ,\ 

s(|a|'  Si  '  S?  ’|B|>  Si 


We  shall  need  the  following  lemmas  : 

Lemma  3.1.  ( Discrete  Sobolev  inequality ) 

Given  e  >  0,  there  exists  a  constant  c  dependent  on  e  aid  n  such  that 

IP^IUoo  <  <IP"«IU2  +  cllulUa  1  <  n> 

||D'«||z.a  <  e\\Dnu\\L3  +  c||u||x,8  l  <  n. 

Lemma  3.2.  (Estimate  ||u"||J. 

For  any  €  >  0  there  exists  a  constant  K  dependent  on  e,  such  that  the  solution  u"  of  (3.1) 
satisfies  : 

IK-'ll2  <  E(«;)2|d2»;i2  +  <ncvn2  +  *(<)  [n«"f  +  urn2] . 

Proof.  From  (3.1),  we  have 

K+1 II2  =  53  {  [*a"f?2«”  +  i(i  -  1  )a1D2v?  +  G?  •  -ia’JD2^  +  i(-*  -  1  )a?D2V?  +  G ) 

i*  l  '  *■ 

where 


Gj  —  BjDuj  +  CjUj  +  fj. 


Expanding  the  above  expression,  we  have 


HOI2  =  tiW'W  +  "  i){^?h2D2^  •  D2^}  +  Ua^h^D^]2} 


+  2 Re  jya]D2u]  +  i(t'  -  l)a]h2D2v^]Gnj 


Using  (3.3),  we  obtain 


-  E  »?l<«"«>2(c”..  - «?)  -  («")!W  -  <5-i)  -  (W+.)2  -  («”)2)(*>“+.  - "")] 

i-i 

=  E  «?-,.;(»?  - «;-)  -  E  WW  - »?-.)  -  E  ”"((‘*”«)2  -  K)2)W«  - 

>■  2  i»l  >“1 

-  -  Ew)2w  -  ^.is  -  E«»”+.>2  -  (°”)2),“(®"+i  -  »;)• 


f* 


hr 


'j  »»■  »■  u*  T7  r >*  ',’  r  r  ’ ."  fy  t  r.«  rr  »-*  r  ’.’v  r  v 
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Hence 

\  f>«[(l  -  0(a?)2^2«"  •  £2^]  + 1  £K)V|i>S"! 


/*1 


i=l 


y-i 


s  -f  D«")!«  -  + 1  I>"+I>!  -  -  <?)l 

>=1  >=1 

+  j£w)>"«-”"i2+w-^1f!] 

>=i 

£  I  X>”+>)2  -  W’KW+i 1  £((«”)!  -  -  'f-il 


>=1 

J-1 


y=i 


So 


s2DK+*)’-K),ii<fii«?«-<fi. 

>«i 

l«?+1||2  <  £(a?)2p2«7|2  +  2  £  l(a,n+1)2  -  (a")2|  K"|  K+1  -  t;"| 

y-i  >=i 


^-1  r 


+  2Re  £  |ta"Z>2«"  +  -  \)anih2D2v^ 

For  C7j ,  using  Lemma  3.1,  we  have 

)|<5"||2  <  3A/2(||Z?«"||2  +  f|«”jj2)  +  3||/"||2 

<  6A/2ei||£>2u"H2  +  (6M2A'(ev)  +  3M2)!|un||2  +  3||/n||2, 
where  «j  >  0.  It  follows  that  if  we  define 


<?>• 


!  1 

Q  =  2ReJ2  *<>"&'<"  +  J(*  -  l)a"h2Z>2t-"  \Gnj, 

i- 1  J 


we  have 


3<6A/j|DV|!  |j<5"j] 


<  3  M 


£2 


3  A/e  2  ■+■ 


18A/3ei 


<2 


|j£>V||2  +  1  [18A/3A (fl)  +  9A/3]  ||un(|2  +  —  ||/n||\ 
£2  «2 


where  €2  is  any  positive  constant.  If  we  define  ei  by 

,2 


<\  - 


86 4M*  ’ 


and  €2  by 


f2  12  M  ' 

then  we  have 

g<i||D2w"||2  +  A'1(0  [||«l8  +  ||/"ll2]. 
where  K\(e)  is  a  constant  dependent  on  e.  On  the  other  hand  using  Lemma  3.1, 

-  («,")2i  ki  k«  -  >-"i  <  m^ikii2 

;=i 

<8M2h[e3||i>2n"||2+A2(e3)||un||2] 

=  8A/2/u3j|Z>V|i2  +  8M2/ili:2(€3)||un||2 

where  €3  is  an  arbitrary  positive  constant.  If  we  take 

6 

€3=  16  APh  ' 

then  we  have 

2  E  IW+i)2  -  («")!l  KI  K«  -  <?l  £  jllcVll2  +  k3MII«”II2  • 

>- 1 

where  A'3(c)  is  a  constant  dependent  on  e.  Combining  the  above  results,  we  finally  have 

K+1II2  <  X>")2 \d2u"\2  +  4D2u"\\*  +  K\ (c) |jun||2  +  ^(cJIIHI8, 

;*=i 

where 

K*(t)  =  K\(c)  +  A’3(e)  , 

which  completes  the  proof. 


Theorem  3.1.  If  (3.4)  holds  and 

min  |a(x)|  >  0  >  0 

then  the  scheme  of  (3.1)  is  stable  if 

1  h2 

”  <  2  max  |a| 

Proof.  Multiply  (3.1)  by  u"  and  we  obtain 

(£«".«")  =  (/",u"). 


(3.12) 


Next  we  estimate  each  term  in  (3.12).  For  the  first  term  we  have 


-  11**11*)  -  yll«r‘llI. 

which  is  similar  to  (2.8).  Applying  Lemma  3.2,  we  obtain,  for  any  e  >  0, 

1  St.  mJ.1  ttO  ■■  Mi.Ox  k 


Re(u^\un)>~(\\u 


The  second  term  is  given  by 


n+lll2-H«"ll2) 


J- 1 


J-i  , 

i=i v 


5>J*)2|Z>2«7|2  +  £|pV||2  +  *<|  |«"| 

;=i 


■  -2«"  +  u"-t  „n\ 


h7 


«?) 


We  first  employ  the  following  expansion  : 


_n  **■)+!  2lly  +  Vj-l  _  JL^  r  _  .  _  _  n\  _  n  /  n  _  n  \1 

ai  ^2  ~  h7  [°/+t/2'l<>+1  ui'  ay— 1/2 u>-i)J 

+  h  [w  -  «"♦■/»><•?«  - «?>  -  W  - 


Then,  using  the  boundary  conditions,  we  have  : 
J-i 


= e  A«*«/tw+.  -  «")*“  -  E - «”)«”+. 

>-i  >-o 


- p  a>+l/2lujf+l  ”  u>|2  » 


J=0 


and 


i  £  p  [(«"  -  «"+i/2)(«"+i  - «")  -  k  -  «;-!/*)(«?  -  «?i 

<  i  E  1 E +|St— 

^  y-i  j*i  /«i 


2  /-I 


I+*lll/"ll 

(3.1 


«*-.)]  - 


<  2[e||D2u"j|2  +  C2||unj|2]  +  JVf2||un|| 
<2€||Z32u"|i2  +  (2c2  +  A/2)||unt|2. 


Hence  we  have 


|*e  J2|  <  2<||0V||2  +  (2c2  +  A/2)|[«"||2. 


The  third  term  is  given  by 


J- 1 

*-£ 

i=i 


_K  J°?l>(v.ti  _  2v"  + •  cpj 


Using  the  expression  : 


(1  -  t)  n  (tf+1  ~  2V"  +  V"_j)  _  (1  -  0  1  r  „  _ 


h 2 


4  fc2 
.  (1  -  0  1 


(a”+lv”+l  -  2alvl  +  al-  lv?-l 


4  fc2 


[(«"  -  a"+iK+i  +  K  - 


we  have 


i,2  *)  \  '  1  r  n  .n  o«n..n  i  „n  -,n  1  -n 

"  2  2^  p  K+Ivi+1  -  2a>  vi  +  a/-itV-iJ  «/ 

y=i 


j-i 


>-i 


and 


£  p  [(<■?  -  <4«K+I  +  <«/  -  «?l 

/“I 


itii  i„" 


A£lJ.a?-a? 


«wi«ji+§£P-t? 


>=i 

<  Mh[e||v”||2  +  -||u 


1"  n  m2i 


e 

ni|2  ,  M^N.,"I|2 


==fc£M|K||2  +  ~||«"|| 


The  rest  of  the  terms  are  easily  bounded  : 


J~l  u -  un 


£  •  «>i  ^  w«i  ikii 

j- 1  in 


M 


<y  [e|pV||2  +  ii:2(e)||«"|i2] 
|(cu,u)|  <  M\\u\\2 


and 


l(/,«)l<;ll/ll!  +  |ll«r  • 


Summing  all  the  above  terms,  we  get 


-  Hu-f) 


2 


J=1  v  7 

[t 


n-fl  *  .  .  ,  avj.  t 

— e  +  2  e  +  heM  +  — - 


IPVII2 


i  -n „  .  ,o  Mh  Affc2  .  1 

+  I  +  2c2  +  A/2  +  ~  +  M  +  - 


Kll2  + 


if  i  ^ 

T*,  +  2 


uni1. 


Of  course,  when  (3.11)  holds,  we  have 


^ki*- 

2  '  J '  4  1  ~  2 


where 


A  = 


h2 


2  max  |a| 


—  />n  >  0. 


As  long  as  e  is  small  enough,  we  can  make 


kn  n  ....  Me  ^  /?2A 

T‘ +  2<  +  +  T  s  V- 


For  such  an  e,  we  have 


-U(K+T  -  Kf)  <  Sf,iiu"ii! + M,nrii!. 


Here  the  constants  Mi,  M2  depend  on  A  and  M.  Using  Lemma  2.1,  we  have  the  proof. 


4.  Concluding  Remarks 

Using  energy  methods,  we  have  established  the  stability  properties  of  the  two  schemes  consid¬ 
ered  in  this  paper.  Since  equation  (1.1)  includes  both  parabolic  and  Schrodinger  type  equations  as 
special  cases,  our  stability  results  provide  a  unified  treatment  for  both  types  of  equations.  Finally, 
the  results  here  agree  with  the  stability  results  obtained  for  the  constant  coefficient  Cauchy  problem 
via  Fourier  analysis  (3,  4]. 
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Appendix 


Appendix  A  (Proof  of  Lemma  2.1) 

We  substitute  (2.5)  into  (2.1)  and  get,  for  the  first  term  : 

«i*+1  =  (<H£+1 

=  l  [r>;+1  - v") + (^n+i  -  *r>r + -  w] 

*  J  [^+°(tj+1  -  t>?)  +  (^n+1  -  *?)((  1  +  a)u”+1  +  otf )] 

S=^+“v?+1  +  <+1-“^+1. 

In  the  following,  for  simplicity,  we  shall  drop  the  index  “n  +  a* .  For  the  second  term  using 
the  expressions 

<t>j+ 1  =  <t>j  +  D(/>jh  +  - D+D~<j>jh 2 


<A>-1  =  4>j  ~  D<j>jh  +  ^D+D-tjh2, 


we  obtain 


D(ADuj)  =  L>(A0(H;) 

=  [Ay+l/2(Vj+l  -  vj)  -  Aj-i/2 (Vj  -  Vj-x)] 

+  2^2  [a/+i/2(^/+i  -  “  •Ai-i/2(^'+i  ~  0y-i)«j-i] 

,  (<Aj+i  _  2<£j  +  ^;-i) , ,  ,  ,  , 

H - 1 - 2^2 - (Aj+l/2vj+l  +  Aj-i/2Vj-i) 

-  %  [^y+l/2(vi+l  -  Vj)  -  Aj-ii2[vj  -  Vy_i)]  +  Aj{Vj+l  -  Vj-i) 

+  [(Ay+1/J  -  Aj)vj+ 1  +  (A,-  -  Ay_x/2)vy_i] 

,  (^/+l  -  2^y  +  <t>j-i) ,  A  ,  ,  x 

+ - 2h? - (Aj+i/2Vj+i  +  Ay_i/2«y_j) 

=  4>jDADvj  +  2AjD4>jDvj  +  Vj+l  +  ^LZ^ZlllVj_l  D<f>j 

,  (Aj+i/2Vj+i  +  Ay_j/jVy_i )  2 

+  - - - p 

For  third  term  using  the  expressions 


we  obtain 


<f>j+ 1  =  <f>j  +  D+Pjh,  <f>j- 1  =  4>j  -  D-<f>jh 


~  B  ■ 

BjDuj  ~2 h^J+lVj+1  ~ 

~  B 

—Bj4>jDvj  +  -£(D+<i>jVj+ j  +  D-tjVj^x)  . 


TTT 


\r\\ 
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In  the  second  and  third  term,  the  first  order  difference  term  for  v  is  (2 AjD<j>j  +  Bj4>j)Dvj.  We 
can  choose  4>j  to  make  this  term  vanish  by  setting 

2AjD<t>j  +  Bj4>j  —  0.  (1) 

Since  |A|  >  ao  >  0  and  A,  B  are  bounded,  it  is  easy  to  see  that  we  can  choose  <j>j  such  that 
0  <  Mo  <  \4>j\  <  Afi,  where  Mo,  Mi  are  constants.  Moreover,  from  (1)  and  \<f>j\  <  Mi,  we  have 

<v 

\D<t>j\  <  constant  in  Qj. 

If  we  take  the  finite  difference  of  (1),  then  because  of  the  assumptions  on  A  ,B,  we  can  obtain 

\D+D~<t>j\  <  constant  in  Qt- 

After  some  manipulation,  we  obtain  (2.6),  where 

l  (An+°  —  An+°  i  i  hn+a 

o"i-°  = 

CTn+a  =  /n»n+o 

i  J 

Kn+a  —  1  b, ^  I  1n  £>  ^n+a^n+cr  ,  h _ fwn+a 

■'+1  4>n+a  h  +  Xi+i/2 +  —D<Pj 

^n+1 

Qn+a  = - k_ 


ff*+a 

pn+a  _  _ 

J  <t>Ta 


Because  of  the  properties  of  A,  B,  C,  F  and  4>,  the  functions  G,  H,  K,  Q  and  F  are  all  bounded. 


Appendix  B 

Lemma  4.1.  (DuhameFs  Principle). 

Assume  un  >  0,  v"  >  0  for  n  =  1,2, _ 


/:nvn,  n  =  0,l. 

•  •  •  1 

n 

un+l  < 

eMr+l  i 
* 

u°  + 

E 

k 

ImO 

and  tn+I 

=  Er.0 

k,. 

Proof. 


.6 


We  have 


Un+1  _|_  j V/fcn)u"  +  knvn 

<(1  +  Affcn)((l  +  A/fc„_i)tin_1  +  fcn-iv"-1]  +  k„vn 


n— 1  n 


<  TT(1  +  Mki)u°  +  y:  JJ  (1  +  Mk,)ktvl  +  k„vn 


1=0  «=(+! 


<  J~J(1  +  Mki)u°  +  ^2  jj(l  +  Mkjkiv1  +  k„v” 


i=o  *=o 


<exp  ln(l  4-  A/fcf)  •  u°  +  y~]  dki 


Since  ln(l  +  Mk/)  <  Mki,  we  obtain  (2). 


Appendix  C 

Lemma  4.2.  Assume  m,  >  0  (*  =  1,  •  •  • ,  n)  and  k  >  0.  The  necessary  and  sufficient  condition  for 
the  polynomial  in  pi 


to  be  nonnegative  definite,  is 


J2  mipf  -  k(j2Pi)* 


E?  ,  -L  ' 

m. 


Proof.  Define  p  =  (pj,  •  •  ■  ,pn)T,  l  =  (1, 1,  •  •  • ,  l)T  and  D  =  diagonol(mi,m-2,-  ■  •  ,m„).  Then  the 
polynomial  in  (3)  can  be  rewritten  as 

pr Dp  -  k{pTl)2  =  pr(D  -  W/r)p 

The  condition  for  the  matrix  Z?  -  kllT  to  be  nonnegative  definite  is  1  -  klTD~ll  >  0  which  is  (4). 
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